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Abstract 

We study the Hamiltonian dynamics of a one-dimensional chain of linearly coupled particles in a 
spatially periodic potential which is subjected to a time-periodic mono-frequency external field. 
The average over time and space of the related force vanishes and hence, the system is effectively 
without bias which excludes any ratchet effect. We pay special attention to the escape of the 
entire chain when initially all of its units are distributed in a potential well. Moreover for an 
escaping chain we explore the possibility of the successive generation of a directed flow based on 
large accelerations. We find that for adiabatic slope-modulations due to the ac-field transient long- 
range transport dynamics arises whose direction is governed by the initial phase of the modulation. 
Most strikingly, that for the driven many particle Hamiltonian system directed collective motion 
is observed provides evidence for the existence of families of transporting invariant tori confining 
orbits in ballistic channels in the high dimensional phase spaces. 

PACS numbers: 05.45.Ac, 05.60.-k, 05.45.Pq, 05.60.Cd 
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I. INTRODUCTION 



Transport phenomena play a fundamental role in many physical systems. In this context 
the so called washboard potential due to its ubiquity and simplicity establishes the proto- 
type of a periodic potential being employed in a number of applications including Josephson 
junctions [1], charge density waves superionic conductors [3|, rotation of drpoles in exter- 



nal fields m , phase- locked loops 
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51] and diffusion of dimers on surfaces 
to quote a few. Applying an external time periodic driving to the washboard potential, 
interesting effects such as phase-locking, hysteresis [isl and stochastic resonance |lj] are 



found. Recent investigations have dealt with the Hamiltonian dynamics of individual par- 
ticles evolving in a washboard potential whose slope is time-periodically varied by a weak 



external monochromatic modulation field 



15|, 
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171]. Remarkably, for the corresponding 



one-and-a-half degree of freedom Hamiltonian system it has been demonstrated that adia- 
batic modulations of the slope of the potential lead to the generation of transient transport 
dynamics related with enormous directed particle flow. Explanation for this behavior has 
been given in terms of the underlying phase space structure of the externally driven one 
degree of freedom system promoting the motion in ballistic channels [l8 ] . Here we consider 
an extension of previous work on individual particle movement and study the case when the 
(many) particles interact due to linear couplings constituting a one-dimensional chain. Mo- 
tion takes place in a phase space for which details of its intricate structures remain elusive 
not at least due to the high dimensionality. Whether in systems with a large number of 
(microscopic) degrees of freedom such (macroscopic) behavior as collective motion leading 
to a directed flow emanates from high dimensional dynamics is not obvious. In systems 
with at least three degrees of freedom irregular dynamics in the form of Arnold diffusion 
arises l^-j22|. This is possible because the KAM tori do not separate the phase space so 
that the various chaotic regions always overlap and in principle orbits can wander along the 
entire stochastic layer permeating the whole phase space. Let us recall that in one-and-a- 
half degree of freedom weakly nonintegrable Hamiltonian systems the motion can stick to 
the hierarchical structure of persisting tori and islands inside stochastic layers being able 
to constitute the above mentioned ballistic channels which support directed transport. The 
question then is whether such trapping regions persist in many degrees of freedom Hamilto- 
nian systems supporting also directed transport. It is the objective of the current work to 
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demonstrate that even in a system of many coupled particles the generation of a directed 
flow going along with collective motion is possible indeed. 



II. THE DRIVEN SYSTEM OF INTERACTING PARTICLES 



We study a one-dimensional chain system consisting of linearly coupled particles with 
Hamiltonian of the following form 



TV r 9 
P 



n=l 



Y + f„(5„) + (7i(5„,() 



n+1 



(1) 



n=l 



wherein p„ and g„ denote the canonically conjugate momenta and positions of the particles 
evolving in the periodic, spatially-symmetric (washboard) potential of unit period, i.e.. 



Mq) = Uoiq + 1) = - cos(27rg)/(27r) . 



(2) 



The external, time- dependent forcing field 



Ui{q,t) = -Fsin(fit + 0o)g 



(3) 



causes time-periodic modulations of the slope of the potential. It has to be stressed that 
there is no bias force involved in the sense that the following average over time and space 
vanishes, i.e. 



T=2n/n 



dq 



dt 



dU{q,t) 
dq 







(4) 



with U{q,t) = Uo{q) + Ui{q,t). The particles interact linearly with coupling strength k. 
Remarkably, as pointed out in prior literature [l^ in the limiting case of uncoupled particles, 
i.e. K = 0, there results an (unexpected) asymmetry of the flux of particles, emanating from 
one potential well, and flowing to the left and right potential wells which indicates the 
existence of directed transport without breaking the reflection symmetry in space and time 
in this system. One reason for the occurrence of phase-dependent directed transport is the 
lowering of the symmetry of the flow in phase space by the ac-field where this asymmetry 



vanishes only for specific values of the initial phase Gq 



15|. In Ref. [16| the authors report 



further on this exceptional situation and show that directed transport is sustained on fairly 
long time scales despite the presence of chaos. In particular it has been demonstrated that for 
sufficiently small forcing frequencies, i7 <^ 1, the width of the arising chaotic layer diverges 
leading to a strong enhancement of the chaotic transport [16]. 

The equations of motion derived from the Hamiltonian in Eqs. ([T]) read as 

+ sin(27rg„) = Fsm{Qt + Qq) + K,[qn+i + Qn-i - 2g„] . (5) 

The analysis in the current paper deals with zero initial phase of the external force term, 
viz. 00 = 0. The influence of the phase Go on the generation of directed motion in the 



periodic potential Uo{q) has been studied in [17|. We stress that averaging over the initial 
phase Bo yields vanishing net current in the periodically driven system. 



III. THE CASE N = 2 

Before we embark on the study of systems comprising of many degrees of freedom let us 
first consider the case of two interacting particles. For a discussion of the dynamics it is 
convenient to introduce the following canonical change of variables induced by the generating 
function: S = \{qi + q2)Px + ~ (l2)Py relating the old and new variables as 

Pl = l{Px + Py), P2 = \{Px~Py). (6) 

Qx = ^(91 + 92), Qj/ = ^(91 -?2)- (7) 

The Hamiltonian expressed in the new variables becomes 

H = i(Pj + _ 1 cos(27rg,) cos{27cQy) 
4 n 

+ 2kQI - 2F sm{nt)Q.^ . (8) 
The corresponding equations of motion are given by 

Qx = — cos(27rQy) sm{27iQx) 

+ Fsm{Qt), (9) 

Qy = - cos(27r(5x) sin(27rQy) 

- 2KQy. (10) 
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Figure 1: Two units: Time evolution of the variable Qx = ((7i + g2)/2 in the driven case oi F = 0.05. 
The remaining parameter values are k = 0.1, Q. = 10^^ and Qq = 0. The initial conditions are 
(7i(0) = 0.337, ^2(0) = 0.347 and pi(0) = ^2(0) = 0. Note the large scale on the ordinate. The 
inset shows bounded chaotic motions of Qx in the undriven case F = 0. 

Evidently, the impact of the external modulation and the linear coupling occurs in sep- 
arate equations. The coupling between the Qx and Qy degree of freedom results from 
parametric modulations of the respective potential force term. Apart from that, equation 
([9]) for the mean value variable Qx is equivalent to that of the system of noninteracting 
particles. Thus, the external modulation may lead to an escape of particles from a potential 
well initiating ongoing rotational motion of the mean value Qx equivalent to the behavior 
observed in the study of noninteracting particles • Distinct to that the amplitude of the 
difference variable Qy is bounded due to the harmonic restoring force associated with the 
K— term. 

In our dynamical studies the initial conditions are chosen such that the units are initially 
contained in a potential well. Furthermore, in the uncoupled case, i.e k = 0, escape from 
the potential well induced by the driving force, as discussed in is excluded. To ensure 
trapping in the driven but uncoupled dynamics the orbits have to lie fairly deep inside the 
separatrix where there is a large island of stability and the dynamics is still regular. 

Likewise, without the weak external ac-field and when the units are coupled the chain is 
supposed to remain trapped in the potential well despite the arising chaotic motion. The 
latter situation is illustrated in the inset in the top panel of Fig. [H The simulation time 



5 




2000 4000 6000 8000 10000 
t 



Figure 2: Illustration of the changes in amplitude of the coordinate Qy = (gi — 52)/2 in the time 
interval < t < 10^. The initial conditions and parameter values are those of the driven case in 
Fig. [TJ Vertical dashed lines mark the instants of time when the amplitude (torus-torus) transitions 
occur. The inset displays the transition from irregular motion prior to escape to quasi-regular 
periodic motion after escape. 



interval corresponds to a hundred times of the period duration of the external modulation 
T = 2ti/VL with VL = 10"^ 

However, when the interacting particles are exerted also to the external ac-force escape 
of the particles from the potential well occurs and results such as those depicted in Fig. [T] 
are obtained. Interestingly, for a slow modulation with frequency Q = 10^^, after escape 
from the potential well taking place at t ~ 800, there occurs long-range directed transport 
of the particles equivalent to the behavior reported in jl^] for the system of noninteracting 
particles. Notice the huge distance covered by the particles. As reflected in the evolution 
of Qy (cf. inset in Fig. [2]) prior to escape the motion is irregular whereas after escape Qy 
oscillates seemingly regularly indicating capture of the motion near an invariant torus (see 
further below). That is, the individual coordinates qi and q2 perform oscillations with a 
mutual phase lag. Moreover, as Fig. [2] reveals the Qy quasi-stable oscillations with a certain 
maximal amplitude are sustained over a specific time interval at the end of which a rather 
abrupt change of the amplitude occurs and the trajectory comes to lie close to another 
torus. On a finer scale one recognizes that the transition from one torus to another one is 
accompanied by a chaotic interlude (see also |23l|). 
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For further illustration of the dynamics in the five dimensional phase space spanned by 
the variables {Px, Py,Qx,Qy,Q = ^t) we utilize Poincare-plots where the cross-section is 
determined as follows: 

^■.= {{Px,Qx,Q)\Qy = 0,Py>0} . (11) 

We represent in Fig. [3]the Poincare-plot in the B — P^— plane corresponding to the dynamics 
of the adiabatically driven system shown in Fig. [H The momentum P^ experiences vast 
alterations in dependence of the value of the phase 0. To be precise, starting from an 
non-inclined potential, i.e. 0o = 0, the force term F{t) = F sm{Qt) produces a negative 
inclination of the potential in the interval < t < tt/Q. It happens that at some instant 
of time tescape < IT / {2Q) the coordinate overcomes the potential barrier and escapes 
from the well. Hence there remains a time interval [tescape, t^/^] during which the particle 
still experiences a force with positive sign which accelerates further the motion in the right 
direction towards increasing values of Q^- For times n/Q < t < 27r/f2 the force acts in the 
opposite direction. In particular for 27r/fi — tescape < t < 27r/fi the momentum evolves with 
its sign reversed compared to the previous acceleration period. With this we can estimate 
the gain in momentum as follows 

APx = \ I + \dTp 



-escape ^ 2,77 /^ — tescape 



^escape ^ 27r/S7 tescape 

X dr [- sin(27rQ^.) cos(27rQj^) + F sin(fir)] . (12) 

For small Q the rapidly oscillating part connected with the first term in the integral averages 
to zero and we find 

F 

APx = 2— COs(fi tescape) ■ (13) 

In general, the smaller Q the higher is the gain in momentum. In principle, for a sufficiently 
small frequency Q the gain can become arbitrarily large (see also With the choice 

Oo = TT for the initial phase of the modulation term the direction of the motion can be 
reversed. In contrast in the undriven case, where the chaotic motion remains trapped in the 
potential well (as depicted in the inset of Fig.[T]), the momentum variable Px accordingly stays 
within the range determined by the extension in momentum direction of the corresponding 
separatrix loop (not shown). 




Figure 3: Poincare-plot presented in the G — P^;— plane. The phase variable G is given mod (27r). 
The parameter values are k = 0.1, F = 0.05, Q = 10"^ and Qq = 0. 

The temporal behavior of the partial energies Ei^2 = 0.5pi 2 + ^0(91,2) is illustrated in 
Fig. m For times t < 800 the energies do not considerably exceed the separatrix level, 
Eg = l/(27r), which is marked by the dashed horizontal line. Escape from the potential 
well takes place at tescape — 800 giving rise to significant increase of the energies reaching a 
maximal value being equivalent to 10^ times the energetic barrier height AE = 1/tt. The 
two partial energies, and likewise the coordinates and momenta, perform oscillations with 
mutual phase lag (see also further above). The inset displays the temporal evolution of 
the energies over one period of the external modulation term where one cannot discern the 
individual trajectories on this scale. 

In terms of phase space structure we recall that in the case of individual particles (i.e. 
K = 0) in the corresponding adiabatically driven system of one-and-a-half degree of freedom 
there arises a broad stochastic layer from the connection of various zones of instability 
due to resonance overlap. Non-contractible KAM tori confine the stochastic layer from 
below and above and form impenetrable barriers for motion in phase space. For the weakly 
nonintegrable system the chaotic sea contains still islands of regular motion. Provided 
these islands possess non-zero winding numbers orbits with initial condition inside such an 
island facilitate transport. Moreover, the motion around these islands is characterized by the 



stickiness to them that can lead to trapping of the trajectory for a long time 



25j. This is 



due the intricate structure of the stochastic layer where close to resonances at the boundary 
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Figure 4: Temporal behavior of the partial energies in the driven system of two coupled units 
around the trapping-detrapping transition. Initial conditions and parameter values as in Fig. [H 
The dashed horizontal line marks the separatrix energy Eg = l/{2ir). The inset illustrates the 
drastic growth of the energies during a period T = 2it/0,. 

between regular and chaotic regions there exists a hierarchy of smaller and smaller islands 
and surrounding cantori. The latter can severely restrict the transport in phase space and 
thus effectively partition the chaotic layer [25]. It seems that the cantori are the less leaky 
the smaller the modulation frequency Q. Hence they form almost impenetrable barriers that 
confine trajectories for a very long but transient period. One should remark that eventually 
this transient period of directed motion terminates because the trajectory escapes through 
one of the holes in the cantori and accesses other regions of the chaotic layer. Therefore the 
motion does not necessarily proceed unidirectionally and unless the trajectory gets captured 



by ballistic channels [iSl] it itinerates within the chaotic layer going along with changes of 
the direction of motion. 

For nonlinear Hamiltonian systems with N >2 degrees of freedom only a few numerical 
results addressing the existence of an enhanced trapping regime are known [26]. It is sup- 
posed that the role played by cantori in driven systems with = 1 is played by families 
of A^— dimensional tori, constituting partial barriers in the 2A^— dimensional phase space, 
where the chaotic trajectory can stick to [26]. On the other hand Arnold diffusion is possi- 
ble and hence in principle a chaotic trajectory, wandering along the entire stochastic layer. 
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can explore the whole phase space [19|,[21j. However, due to stickiness to higher dimen- 
sional invariant tori Arnold diffusion can be suppressed so that certain stochastic regions 
are distinguished in which the trajectories become trapped for longer times 3]. 

In fact, our findings imply that in the driven two-particle system the motion takes place in 



ballistic channels 



ISj associated with stickiness to two-dimensional invariant tori (cf. Fig. [3]). 



Notably, in the two particle system the mean value Qx evolves in the same manner as the 
single coordinate in the individual particle counterpart, viz. it exhibits effective growth. We 
performed numerical studies for a large number of initial conditions contained in a potential 
well and found behavior equivalent to that illustrated in Figs. [THU The slight alterations in 
the slope of the graph of in Fig. [T] can be explained with the fact that after a while the 
trajectory can get released from a torus region to be captured in another one where it again 
sticks to for a certain time. As noted above this release and capture behavior is manifested 
in the changes of the maximal amplitudes of the oscillations of the variable Qy (see Figs. [1] 
and El). 

Alternatively to the geometrical phase space interpretation of the occurrence of directed 
motion supported by ballistic channels it is worthwhile to mention that from an analytical 
point of view the so called gating mechanism can be responsible for the initiation of the 
motion [2?! . This mechanism is based on an interplay between an ac driving force and an ac 
parametric force causing dynamical symmetry breaking that may lead to directed motion. 
In fact Eq. ([9]) for the mean value variable Qx exhibits both features, i.e. it contains an ac 
driving term as well as the ac parametric force suggesting relevance of the gating mechanism. 



IV. EXTENDED COUPLED OSCILLATOR SYSTEMS 

Finally we consider extended coupled oscillator chains where we aim to demonstrate that 
it is also possible to generate directed motion with adiabatic periodic modulations of the 
slope of the spatially periodic potential. 

Concerning the initial conditions we proceed as follows: Initially an amount of energy 
En = 0.5pn + Uo{q„) < AE is applied per unit such that the whole chain is elongated 
homogeneously along a fixed position go near the bottom of the well. Then, the position 
and/or momenta of all units are randomized. The random position values are chosen from a 
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Figure 5: Typical spatio-temporal evolution of the coordinates Qnit) for a chain of 100 coupled 
particles. Initially the coordinates are uniformly distributed within the interval |^n(0) — qo\ < Aq 
with an average % = —0.35 and width Aq = 0.01 and zero momenta, i.e. Pn(0) = Po = Ap = 0. 
The coupling strength is k = 0.3 and the remaining parameter values are F = 0.05, Q = 10^'^ and 
Oo = 0. Top (Bottom) panel: = 10"^ {Q = 10"^). 

bounded interval |q'n(0) —qo\^ and, likewise, the random initial momenta, |pn(0) — Po| < 
Ap. The whole chain is thus initialized close to an almost homogeneous state, but yet 
sufficiently displaced {Aq ^ 0) in order to generate nonvanishing interactions, entailing the 
exchange of energy among the coupled units. 

First we note that in the case without external modulation of the slope there occurs the 
formation of a pattern of localized states due to modulational instability (not shown here). 
Due the irregular dynamics it happens that occasionally a unit overcomes the potential 
barrier but no coordinated motion of the chain results. 

Remarkably, applying the adiabatic modulation the entire chain not only escapes from the 
potential well but manages also to travel freely and unidirectionally over a giant distance 
as seen from the spatio-temporal evolution of qn{t) in Fig. [51 In comparison for faster 
modulations no directed motion of the chain is obtained. The chain consists of 100 coupled 
oscillators and open boundary conditions are imposed. The profile of the chain continuously 
undergoes changes with ensuing deviations from a flat state. Nevertheless the intriguing 
feature of transients of extremely long-range transport of the chain is provided by collective 
motion which is also reflected in the temporal behavior of the mean value of the coordinate < 
q >=< jj: J2n=i In > shown in Fig. [6l Ensemble averages, denoted by < > were performed 
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Figure 6: Temporal evolution of the mean value < q > for systems driven with different driving 
frequencies as indicated in the plot. Ensemble averages, denoted by < > were performed over 100 
realizations of trajectories whose initial condition lie in the range given in Fig. El The parameter 
values are those given in Fig. \5\ 



over 100 realizations of trajectories whose initial condition lie in the range given in the 
caption of Fig. [51 Note that while for driving with Q < 10^^ on average no directed motion 
results for sufficiently slow modulations (here illustrated for Q = 10~^) a large directed 
current is observed. Conclusively, concerning transport the many particle system exhibits 
features equivalent to that observed in the system of individual particles IGj, 17|. 

It is interesting to relate our findings to those obtained for the continuum case of Eq. (I5l), 
viz. a perturbed sine-Gordon equation that is driven by a mono-frequency periodic external 
field destroying integrability 28|. There also the emergence of a directed current for driven 
solitons results for suitable parameter values of the external field. In comparison we note 
that in the ID continuum case a localized wave-like excitation in the form of a soliton moves 
along the spatial coordinate whereas in our present study of the discrete counterpart the 
entire chain is coherently transported over the wells of the periodic potential landscape. 
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V. SUMMARY 



To summarize, we found that in the time-dependent driven Hamiltonian dynamics of a 
chain of hnearly interacting particles evolving in a symmetric, spatially periodic potential di- 
rected transport can be induced. That is, when the slope of the potential is time-periodically 
varied due to a slowly oscillating external modulation field the transport proceeds unidirec- 
tionally in a long-lasting transient period such that the chain covers huge distances. The 
behavior of the mean value of the coordinates of extended chains is strikingly similar to that 
found in the system of noninteracting particles. Thus, the collective directed motion of the 
numerous microscopic degrees of freedom is manifested at a collective level in the evolution of 
the mean value of q. We recall that for the system of individual particles transporting island 
structures corresponding to ballistic channels present in mixed phase space are responsible 
for transients of directed long-range transport. Our findings for the many particle system 
imply that ballistic channels exist in a high dimensional phase space too and support the 
collective motion of the coupled oscillators of a chain. For systems involving many degrees 
of freedom the issue of stickiness of invariant tori in high-dimensional phase space needs to 
be explored in more detail (an analysis which goes beyond the scope of the current article). 



[1] A. Barone and G. Paterno, Physics and Applications of the Josephson Effect, (Wiley, New 
York, 1982). 

[2] G. Gruner, A. Zawadowski, and P.M. Chaikin, Phys. Rev. Lett. 46, 511 (1981). 

[3] P. Fulde, L. Pietronero, W.R. Schneider, and S. Strassler, Phys. Rev. Lett. 35, 1776 (1975). 

[4] D. Reguera, J.M. Rubi, and A. Perez-Madrid, Phys. Rev. E 62, 5313 (2002). 

[5] C.W. Lindsey, Synchronization Systems in Communication and Control (Prentice-Hall, En- 

glewood Cliffs, New Jersey, 1972). 
[6] J.W.M. Frenken and J.F. Van Der Veen, Phys. Rev. Lett. 54, 34 (1985). 
[7] O.M. Braun, Surf. Sci. 230, 262 (1990). 
[8] E. Pijper and A. Fasolino, Phys. Rev. B 72, 165328 (2005). 
[9] O.M. Braun, R. Ferrando, and G.E. Tommei, Phys. Rev. E 68, 051101 (2003). 
[10] E. Heinsalu, M. Patriarca, and F. Marchesoni, larXiv:0708!2858V l 



13 



[11] M. Partriarca and P. Szelestey, Act. Phys. Pol. 36, 1745 (2005). 
[12] C. Fusco and A. Fasolino, Thin Solid Films 428, 34 (2003). 

[13] M. Borromeo, G. Constantini, and F. Marchesoni, Phys. Rev. Lett. 82, 2820 (1999). 
[14] L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). 
[15] O.M. Yevtushenko, S. Flach and K. Richter, Phys. Rev. E 61, 7215) (2000). 
[16] S.M. Soskin, O.M. Yevtushenko, and R. Mannella, Phys. Rev. Lett. 95, 224101 (2005). 
[17] D. Hennig, L. Schimansky-Geier, and P. Hanggi, Eur. Phys. J. B 62, 493 (2008); Europhysics 
News 39, 21 (2008). 

[18] S. Dcnisov and S. Flach, Phys. Rev. E 64, 056236 (2001); S. Dcnisov, J. Klaftcr, M. Urbakh, 
and S. Flach, Physica D 170, 131 (2002); S. Denisov, J. Klafter, and M. Urbakh, Phys. Rev. 
E 66, 046217 (2002). 

[19] V.I. Arnold, Sov. Math. Dokl. 2, 245 (1961); ibid 5, 581 (1964). 

[20] V.I. Arnold Mathematical Methods of Classical Mechanics (Springer, New York, 1978). 
[21] B.V. Chirikov, Phys. Rep. 52, 263 (1979). 

[22] A.J. Lichtenberg and M.A. Lieberman, Regular and Stochastic Motion (Springer- Verlag, New 

York, Heidelberg, Berlin, 1983). 
[23] G. Zumofen and J. Klafter, Europhys. Lett. 25, 565 (1994). 

[24] M. Zaslavsky, Chaos in Dynamical Systems (Harwood, New York, 1985), Physics of Chaos in 

Hamiltonian Systems (Imperial College Press, London, 1998). 
[25] R.S. MacKay, J.D. Meiss, and I.C. Percival, Physica D 13, 55 (1984); J.D. Meiss and E. Ott, 

Physica D 20, 387 (1986). 
[26] H. Kantz and P. Grassberger, Phys. Lett. A 123, 437 (1987); M. Ding, T. Bountis, and E. 

Ott, iUd, 395 (1990); H. Kantz and E.G. Altmann, Europhys. Lett. 78, 10008 (207). 
[27] E. Zamora-Sillero, N.R. Quintero, and F. Mertens, Phys. Rev. E 76, 0666601 (2007). 
[28] O.H. Olsen and M.R. Samuelsen, Phys. Rev. B 28, 210 (1983). 



14 



